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SUMMARY 


Acid-base  concentration  cells  using  room  temperature  ionic  liquids  have  been 
constructed,  discharged,  and  analyzed  in  previous  work  in  the  Frank  J.  Seiler 
Research  Laboratory.  The  ionic  liquid  electrolytes  in  these  cells  are  formed 
by  mixing  l-methyl-3-ethylimidazolium  chloride  (MCI)  with  AlCl^  at  various 
mole  fractions,  N,  of  AlCly  The  proposed  electrode  processes  in  the  cell 
are  oxidation  of  the  aluminum  anode  in  basic  melts,  N  <  0.5,  and  reduction 
of  A^Cl^  in  acidic  melts,  N  >  0.5,  to  A1  at  the  cathode.  The 
theoretical  specific  energies  of  such  cells  depends  on  the  values  of  N  for  the 
acidic  and  basic  melts  and  the  transport  numbers  of  the  M+,  AlCl^  ,  and 
Al^Cl^  ions.  This  report  shows  how  the  specific  energy  depends  on 
these  parameters  and  how  theoretical  specific  energies  were  calculated. 
Theoretical  results  are  also  shown  to  agree  with  previously  obtained 
experimental  results. 


PREFACE 


This  is  a  report  on  work  done  in  Work  Unit  2303-F2-11,  Room  Temperature 
Molten  Salt  Battery  Research.  The  report  is  intended  to  show  how  calculations 
of  specific  energy  made  in  conjunction  with  experimental  work  done  in  this 
work  unit  were  actually  carried  out.  Results  of  the  experimental  work  are 
being  reported  in  the  open  literature,  unburdened  by  the  relatively  laborious 
step-by-step  presentation  of  the  calculation.  I  wish  to  thank  Capt  Auborn, 
USNR,  and  Bell  Research  Labs;  Lt  Col  Fannin,  Maj  Williams,  Dr  Wilkes,  and  Dr 
Donahue  of  the  Frank  J.  Seiler  Research  Lab;  and  Dr  Lowell  A.  King  for  many 
helpful  discussions  regarding  the  validity  and  usefulness  of  these 
calculations.  I  also  appreciate  the  invaluable  assistance  of  Lt  Col  Fannin  in 
the  programming  of  the  calculations.  I  especially  thank  Judy  Cross  for  her 
cheerful,  persevering,  and  expert  secretarial  assistance  in  the  preparation  of 
this  report. 
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I.  Introduction: 


Researchers  in  this  Laboratory  have  discovered  a  chloroaluminate  salt 
which  is  liquid  at  or  well  below  room  temperature  (1).  This  salt  is  formed  by 
mixing  l-methyl-3-ethylimidazolium  chloride  (MCI)  with  AlCl^,  both  solids  at 
room  temperature,  as  shown  below: 


F~\  ci- 
ch^  c2ti5 

/7\ 

+  Alcl3  a1c1a  *  K>>1- 

(1) 

(MCI) 

(M+) 

When  AlCl^  is  added  in  less  than  stoichiometric  amount,  i.e.,  the  mole 
fraction  of  AlCLj  in  the  melt,  N,  is  less  than  0.5,  the  anions  in  the  melt 
are  AlCl^  and  Cl  .  Such  mixtures  are  basic  since  Cl  is  a  Lewis 
base.  When  N  >  0.5,  the  process. 


aici3  +  aici4  — ►  ai2ci?  ,  K»l, 


(2) 


occurs,  and  the  anions  are  AlCl^  and  AljCl^  .  Such  melts  are  acidic 
since  A^Cly  is  a  strong  Lewis  acid.  The  anion  mole  fraction  of  each 
ion  as  a  function  of  N  is  shown  in  Fig  1.  At  mole  fractions  greater  than 
0.50,  the  reaction. 


ai2ci?  +  aici3  - ►  a13C110’ 


(3) 


may  also  be  occurring;  but  we  shall  ignore  it  in  the  present  discussion. 

The  acid-base  properties  of  these  melts  are  best  understood  in  terms  of 


the  equilibrium. 


A12C1?  +  Cl  2A1C14  ,  K»l.  (4) 

Since  K  for  this  acid-base  neutralization  reaction  is  large,  acidic  and  basic 
melts  can  be  combined  in  concentration  cells,  as  has  been  demonstrated  (2). 

In  these  cells,  a  schematic  diagram  of  which  is  shown  in  Fig  2,  measured 
specific  energies  were  shown  to  approach  calculated  specific  energies.  The 
purpose  of  this  report  is  to  show  how  the  specific  energy  of  such  a 
concentration  cell  can  be  calculated,  given  N's  of  the  acidic  and  basic  melts 
and  transport  numbers,  t±  for  the  M+  and  Al2Cl7~  ions.  It  can 
be  shown  that  the  transport  number  for  AlCl^”  will  not  affect  the 
specific  energy  for  a  given  t_( Al 2C17~) •  This  can  be  seen 
intuitively  from  Fig  2  which  shows  that  of  the  three  transporting  ions,  only 
Al2Cly  is  involved  in  the  electrochemical  reaction.  Therefore,  in 
this  approach,  t+(M+)  is  treated  as  effectively  the  sum  of  t+(M+)  and 
t_(AlCl4  ),  so  it  is  not  the  same  as  the  true  transport  number  for  M+ 
obtained  experimentally.  Since  Cl”  flow  out  of  the  anolyte  would  oppose  the 
current  and  therefore  should  not  occur,  t_  (Cl”)  is  set  at  zero. 

II.  Theory: 

Approach.  The  anode  and  cathode  reactions  shown  in  Fig  2  do  not  take  into 
account  the  mechanism  for  ion  transport.  Assuming,  based  on  the  discussion 
above,  that  only  M  and  A12C17  transport  charge,  the  unbalanced 
anode  and  cathode  reactions  are: 


MOLE  FRACTION  AICI3 

Anion  mole  fractions  in  chloroaluminate  melts  of  various  AICI3  mole 
fraction. 
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g  2.  Schematic  diagram  of  concentration  cell. 


(5) 


ANODE:  (MCI  +  AlCl.^)  +  A1  +  A12C17~  — >  MAlCl4  +  M+  +  e 

CATHODE:  (MCI  +  A1C13)  +  M+  +  e“  — *  Al  +  MA1C14  +  A12C17“  (6) 

These  reactions  differ  from  those  in  Fig  2  in  that  1)  M  and  Al2Cl7 

are  included  so  that  transport  of  charge  by  these  species  can  be  accounted  for 

and  2)  the  actual  reacting  species  in  the  starting  electrolytes  (MCI  + 

A1C13)  are  not  specified,  thus  allowing  the  mass  balance  to  be  carried  out 
for  various  to-be-specified  mole  fractions  of  A1C13 .  In  each  equation,  the 
(MCI  +  A1C13>  represents  the  starting  melt  components  with  mole  fractions  to 
be  specified.  MA1C14 ,  or  N  =  0.5  melt,  is  on  the  product  side  since  in 
theory  the  cell  can  be  discharged  until  both  acid  and  basic  melts  are  neutral 
if  stoichiometric  amounts  of  these  melts  are  used  as  the  catholyte  and 
anolyte.  The  approach  then  is  to  specify  the  starting  mole  fractions  of 
anolyte  (NA)  and  catholyte  (NC)  and  the  transport  numbers  for  M+  and  Al2Cl7 
(subsequently  designated  t+  and  t_7)  and  then  balance  the  equations  for  a 
given  amount  of  charge  passed.  From  the  coefficients  of  the  MCI  and  A1C13 
in  the  acid  and  basic  melts  and  the  Al  reacted,  we  get  the  masses  of  anolyte 
(MALY)  and  catholyte  (MCLY)  and  Al  transferred  (MAL)  required  to  pass  the 
given  amount  of  charge.  Then  specific  energy  is  calculated  from  equation  (7) 

(V)  (Q/96,487)eq  (26.8  Ahr/eq) 

SP  EN  (Whr/kg)  =  -  (7) 

(MALY  +  MCLY  +  MAL)  kg 

where  V  is  the  load  voltage,  Q  is  the  charge  passed,  and  26.8  Ahr/eq  is  a 
conversion  factor.  All  the  terms  in  eqn  (7)  can  be  measured  experimentally 
(2).  Thus  the  specific  energies  so  calculated  can  be  compared  to  experimental 
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values  to  test  the  model  represented  by  Fig  2  and  equations  (5)  and  (6). 

2.  Determination  of  Coefficients  in  Balanced  Cell  Reactions. 

a)  We  start  by  specifying  that  in  the  theoretical  cell,  one  equivalent  of 
charge  is  passed,  i.e.,  the  coefficient  of  e  is  1. 

b)  Next  we  assign  a  value  for  t_7  which  fixes  t+  at  l-t_7  (recalling 

that  t  is  effectively  the  sum  of  t  and  t  ,  and  that  1  =  t.  +  t  ,  + 

+  +  -4  +  -4 

t_7). 

c)  When  the  t_7  A^Cl^  is  transported  across  the  separator  into  the 
basic  melt  which  contains  MCI,  the  following  reaction  occurs. 


& 

ft 

& 


t_7  A12C1?  +  t_7  MCI 


t  MA1C1.  +  t  A1C1. 
-7  4-74 


(8) 


in  accordance  with  equation  (4).  Since  it  doesn't  matter  when  in  the  course 

of  the  cell  discharge  this  occurs,  we  can,  for  the  purpose  of  simplifying  the 

mass  balance,  replace  the  reactant  Al2Cl7  in  the  ANODE  reaction  by  the 

t  _  A1C1.  and  t  MA1C1.  formed. 

-74  -7  4 

d)  Ail  the  available  evidence  (2)  demonstrates  that  1/3  mole  Al  is  consumed 
at  the  anode  and  deposited  on  the  cathode  for  each  mole  e  passed.  Thus  the 
coefficient  of  Al  in  both  reactions  is  set  at  1/3. 

e)  At  this  point,  the  ANODE  reaction  can  be  written. 


ANODE:  (MCI  +  A1C1,  )  +  1/3  Al  +  t  ,  A1C1.  +  t  ,  MA1C1.  ,  , 

3  _  -7 _ 4  -7 _ 4  ( r; 

MA1C1.  ,  .  +  t  M+  +  e" 

4  (p)  _+ _  _ 


(9) 


where  the  underlined  species  have  their  coefficients  specified.  The  (p)  and 
(r)  subscripts  distinguish  the  reactant  MAlCl^  formed  by  equation  (8)  from 
the  product  MAICI^  of  the  electrochemical  reaction.  Now,  balance  of  Al  in 


equation  (9)  requires  that 


5 


(10) 


C(A1C1^ )  =  C(MAlCl4(p))  -  1/3  -  2t 

and  balance  of  Cl  requires  that 

4  C(MAlCl4(p))  -  C(MC1)  +  3  C(AlCl3)  +  8  t_? ,  (11) 

where  the  C's  are  the  coefficients  of  the  designated  species  in  the  ANODE 
reaction. 

f)  The  MCI  in  equation  (9)  is  what  remains  of  the  MCI  initially  in  the  melt, 
MCl^,  after  reaction  with  the  transported  Al^l^  (equation  (8)). 

Thus  we  can  write, 

C(MCli)  ■  C(MC1)  +  t_?.  (12) 

Now  C(MCL )  is  related  to  NA,  the  initial  mole  fraction  of  AlCl^  in  the 
anolyte  melt,  by 

C(MC1. )  ■  Il-Z  NA>  C(A1C1_)  (13) 

NA 

Combining  equations  (12)  and  (13)  gives 
C(MC1)  ■  il  ~  NA>  C(A1C1.)  -  t  , 

NA  3  _7  (14) 

Substitution  of  equation  (14)  into  (11)  gives 


4  C(MAlCl4(p))  =(1  ~  NA)  C(A1C13)  -t_7  +  3C(A1C13)  +  8t_7 


or  simplifying, 


4C(MA1C1.  ,  . )  =  3  +LI_^NA)  C(A1C1,)  +  7t  7  (] 

4(p)  NA  3  ~7 

If  we  then  substitute  the  expression  for  C(A1C1.  )  in  eqn  (10)  into  (15),  we 


4C(MA1C1.  .  .  )  =  3  +(1_I_NA)  C(MA1C1.  ,  .)  -  1/3  -  2t  +  7t  _ 

4(p)  [  NA  J[  -7J  *7, 


which  can  be  rearranged  to  give 


,  ,  -  |~3+(1  "NA)1  (1/3 

C(MA1C1  ,  )  = - 1 M - J - 

4(p)  4  -  3+11 


(1/3  +  2t_?  )  +  7t_ 
4-  (1  "  NA)] 


This  is  the  equation  used  to  calculate  the  coefficient  of  MAlCl^^j  in  the 

ANODE  reaction  for  specified  values  of  t_7  and  NA. 

g)  Now  that  we  have  C(MAlCl^^j),  we  get,  from  balance  of  M, 


C(MCli)  *  C(MAlCl4(p))  +  t+. 


and  with  C(MCl^)  determined,  we  also  have  C(MCl)  from  equation  (12)  and 


C(AICI-)  from  eqn  (13). 


At  this  point,  all  coefficients  in  the  ANODE  reaction,  eqn  (9),  have 
been  calculated,  including  the  initial  number  of  moles  of  MCI,  C(MCl^), 
required  to  produce  1  equivalent  of  electrons  by  stoichiometric  reaction  to 


the  neutral  melt. 

h)  By  the  same  arguments  as  given  in  steps  (a)  through  (e),  we  can  write  the 
cathode  reaction. 


CATHODE:  (MCI  +  A1C13>  +  t+M  +  le  -•>  1/3  A1  +  MAlCl^  +  t_?  Al2Cl?  (18 


where  the  coefficients  of  the  underlined  species  are  known.  Thus,  we  need  to 
determine  the  CATHODE  reaction  coefficients,  B(MCl),  BCAICI^),  and  B(MAlCl^) 
for  a  specified  mole  fraction  of  the  catholyte,  NC. 

i)  We  start  by  noting  that  the  amount  of  A^Cl^  which  must  be  present 
in  this  melt  to  produce  1  mole  e  by  the  CATHODE  reaction  shown  in  Fig  2  and 
transport  t_^  moles  of  charge  is  (4/3  +  t_^).  Then  we  ask  how  much 
AlCl^  and  MCI  must  be  present  in  the  melt  of  mole  fraction,  NC,  to  contain 
this  amount  of  A^Cl^  .  Since  the  reactions, 


MCI  +  A1C1,  +  MA1C1,  , 
3  4’ 


maici4  +  aici3  +  mai2ci?, 


both  have  K»l,  we  can  write 


mol  (AljjClj  )  *  mol  (A1C13)  -  mol  (MCI),  or 


(4/3  ♦  t„7)  -  NC  n  -  (1  -  NC)  n  -  (2NC-l)n 


where  n  is  the  total  moles  (MCL  +  AlCLj)  of  catholyte.  Thus,  since 
B(A1C13)  *  NC  n,  we  can  write,  using  eqn  (21), 


M* 

M* 

Pv 

•v; 


p 

i 


B(A1C13)  -  NC 


Similarly, 


4/3  ♦  t_7 


2NC  -  1 


(22) 


I  B(MC1)  -  (1  -  NC) 

4/3  +  t-7 

4 

2NC  -  1 

(23) 


j)  Balance  of  M  in  eqn  (18)  requires  that 
B  (MAICI^)  -  B(MC1)  +  t+ 


(24) 


Thus  all  coefficients  in  the  CATHODE  reaction  to  produce  1  mol  e  and 
neutral  melt  are  also  determined  for  given  values  of  NC  and  t+. 
k)  If  we  specify  NA,  NC,  and  the  load  voltage,  V,  we  can  now  calculate  the 
specific  energy  for  a  range  of  values  of  t  by  using  eqn  (7)  in  the  form, 


SP  EN  *  £(26.8)  x  (V)]/ Jo.  14664  (B(MC1)  +  C(MCl.)) 


+  0 . 13334 (B( A1C1  _)  +  C( A1C1,))  ♦  P.t.02698 
3  3  3 


(25) 


where  the  coefficients  in  the  denominator  are  multiplied  by  the  appropriate 
molecular  weights  in  kg.  This  equation  appears  in  the  program  shown  in  Fig  3 
as  line  1240.  This  is  a  program  which  uses  our  Laboratory's  Program  for 
Handling  Data  (PHD)  described  in  Ref  (3).  It  performs  the  calculation  of 
coefficients  as  described  above. 

1)  Since  it  is  not  practical  to  construct  cells  with  exactly  stoichiometric 
amounts  of  anolyte  and  catholyte,  we  have  chosen  to  make  the  anolyte  the 
limiting  reagent.  To  take  this  into  account  in  theoretical  calculations  of 
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SPEHPHO  20-OCT-1983  14IS6 

990  INPUT  ‘HOLE  FRACTION  OF  ANOLVTE,  NA'jNA 

991  INPUT  '  HOLE  FRACTION  OF  CATHOLVTE,  NC'jNC 

992  INPUT  ‘BATTERY  LOAD  VOLTAGE,  V,V 

993  INPUT  ‘ELECTROLYTE  NASS  RATIO.  HASS  CLY/HASS  ALY,  ELECRAT * j ELECRAT 

995  INPUT  ‘FILE  NANE‘|FILE0 

996  OPEN  FILES** .DAT*  FOR  OUTPUT  AS  FILE  02 
990  PRINT  S2, *?, *|FILES 

999  PRINT  •2,‘NA,.0‘sPRINT  S2. ‘NC, ,0‘sPRINT  #2, ‘T?,,0‘\PRINT  #2.‘TH,.0‘ 

1000  PRINT  02, ‘V, UOLTS, 0‘sPRINT  #2, ‘SPEN.UHR/KG.O* 

1001  PRINT  02 , * SPENXCLY , UHR/KG , 0  * 

1002  FOR  TH-.5  TO  1.00002  STEP  .05 

1003  T7-1-TH 

1010  PRINT  ‘  TRANSPORT  NUHBER  FOR  AL2CL7-.T7* jT7 

1020  PRINT  'TRANSPORT  NUHBER  FOR  IHIDA20LIUN  ION,  TH* jTH 

1080  CH-THNC4»T7\CHACR-T7\CAL-l/3 

1090  CHACP-(-(3*(l-NA)/NA)«(CAL*2*T7)*7tT7)/(4-(3*(l-NA)/NA)) 

U00  cncLi-cn*CHACP 

1110  CHCL-CHCLI-T7 

1120  CALCL3a(NA/( 1-NA ) XCHCLI 

1150  BAL-1/3SB7«T7\BH-TH 

1160  BALCL3-NC* ( ( T7+4/3  )/<  20NC-1 > ) 

1170  BHCL- ( 1-NC )*< (T7+4/3 )/<2«NC-l > ) 

1180  BHAC'BHCL+BH 

1190  ELECRATS-(BALCL3*133.34*BflCL*146.64>/<CHCLI»146.64+CAlCL3*133.34> 

1195  EXCLYRAT-ELECRAT-ELECRATS 

1200  EXCLY-EXCLYRAT*<CHCLI* . 14664+CALCL3*. 13334 ) 

1220  PRINT\PRINT\PRINT 

1230  PRINT  'GIVEN  THE  ABOVE  PARAHETERS* ‘SPRINT 

1240  SPEN-26.8*V/<0.14664*<BHCL*CHCLI>*0.13334*<BALCL3*CALCL3>*0.02698«(CAL>> 

1243  IF  EXCLY<0  THEN  SPENXCLY-Os  GOTO  1250 
1245  SPENXCLV-26 . 8«V/ ( (26 . 8IV/SPEN ) *EXCLY ) 

1250  PRINT  ‘BALANCED  ANODE  REACTION  IS‘\PRINT 

1260  PRINT  USING  •< tt.tttt  HCL  *  tt.tttt  ALCL3)  *  MU  AL  *  0.0000  ALCL4-  *  Mt 
t#  HALCL4*.  CHCL,CALCL3,CAL,C4,CHACR 

1265  PRINTS  PRINT  USING  *— >  tt.tttt  HALCL4  ♦  0.0000  N+  ♦  1.0000  E-‘.CHACP,CH 
1270  PRINTsPRINTsPRINT  ‘BALANCED  CATHODE  REACTION  IS'sPRINT 

1280  PRINT  USING  ‘(tS.Sttt  HCL  *  U.lttt  ALCL3)  *  t.UU  H*  *  1.0000  E-*,BHCL,BAL 
CL3.BH 

1290  PRINTsPRINT  USING  ' - >  MU  AL  *  tt.tttt  HALCL4  *  # . 0BB0AL2CL7- * , BAL, BHAC, 

B7 

1292  PRINTsPRINTsPRINT  USING  ‘THE  STOICHIOHETRIC  RATIO  OF  HCLV/HALY,  ELECRATS.  IS 
00. 000 ‘.ELECRATS 

1294  PRINT  SPRINT  USING  ‘THE  SPECIFIC  ENERGY  IS  ttt.ttt  UHR/KG’, SPEN 

1296  IF  EXCLV<0  THEN  PRINTsPRINT  USING  'THE  STOICHIOHETRIC  RATIO  EXCEEDS  THE  SPEC 

IFIED  HCLY/HALY  RATIO,  t. tf .ELECRAT 

1300  IF  EXCLY<0  GOTO  1302 

1301  PRINTsPRINT  USING  ‘SPEC  ENERGY  INCLUDING  EXCESS  CATHOLVTE  IS  00.00  UHR/KG*, S 
PENXCLV 

1302  PRINT  02, '2* , ';NA; ' , * jNCj  *, ‘ jT?j  ‘ , * }TH; * , * jVj ‘ ,  * jSPEN; * , ‘ ;SPENXCLV 
1305  NEXT  TH 

1307  CLOSE  02 
1310  END 

Ready 


SPENPHD  program  for  calculating  concentration  cell  specific  energy 
for  a  range  of  transport  numbers. 


the  specific  energy,  we  can  specify  an  actual  ratio  of  mass  catholyte  to  mass 
anolyte  (ELECRAT  in  the  program).  The  corresponding  stoichiometric  ratio, 
ELECRATS,  is  computed  directly  from  the  coefficients  of  the  melt  components  in 
the  anolyte  and  catholyte  obtained  as  described  above.  The  mass  of  catholyte 
used  in  the  stoichiometric  cell  is  simply  ELECRATS  x  mass(anolyte) ,  and  the 
mass  of  catholyte  used  in  the  cell  with  excess  catholyte  is  ELECRAT  x  mass 
(anolyte).  Thus,  the  difference  in  these  two  masses,  EXCLY  (calculated  in 
line  1200)  is  simply  added  to  the  denominator  in  eqn  (25)  to  get  SPENXCLY,  the 
specific  energy  for  a  given  mass  ratio  of  catholyte  to  anolyte,  in  which  there 
is  excess  catholyte.  If  ELECRATS  is  larger  than  ELECRAT,  so  that  EXCLY  is 
negative,  SPENXCLY  is  not  calculated  but  set  equal  to  zero. 

III.  Results: 

Results  for  a  cell  with  t_^,  t+,  NA,  NC,  V,  and  ELECRAT  all  specified  can 
be  obtained  using  the  SPENCALC  program  shown  in  Fig  4.  This  program  is 
virtually  identical  to  SPENPHD  except  only  one  set  of  transport  numbers  is 
used  and  the  results  are  not  transferred  into  a  dataset  in  the  PHD  files. 

The  results  for  a  given  run  of  SPENPHD  which  is  representative  of  actual 
concentration  cells  discharged  in  this  Laboratory  (2),  is  shown  in  Fig  5. 

Results  of  runs  done  with  SPENPHD  are  plotted  in  Figures  6  and  7.  In 
Fig  6,  the  catholyte  in  each  run  is  NC  *  0.67,  and  NA  is  varied.  The  specific 
energies  including  excess  catholyte  (SPENXCLY)  are  shown  as  dashed  lines. 

Note  that  there  is  no  curve  of  SPENXCLY  for  the  NA  =  0.03  case  because 
ELECRATS  for  that  case  exceeds  the  specified  ELECRAT  of  1.5. 

In  Fig  7,  the  runs  have  the  same  NA,  set  at  0.33,  and  NC  is  varied. 
Note  that  all  SPENXCLY'S  are  the  same  since  the  denominator  in  the  calculation 
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SPENCALC  29-0CT-1983  14 l S3 

»  1919  INPUT  *  TRANSPORT  NUNIER  FOR  AL2CL7-,T7'jT7 

1020  INPUT  'TRANSPORT  NUNIER  FOR  INIDAZOLIUN  ION,  TH'jTN 
1030  INPUT  'HOLE  FRACTION  OF  ANOLVTE,  NA'jNA 
1040  INPUT  'HOLE  FRACTION  OF  CATHOLVTE,  NC'jNC 
IASI  INPUT  'IATTERV  LOAD  UOLTACE,  U‘|U 

lOSt  INPUT  * EXPERIMENTAL  RATIO.  NASS  CLV/NASS  ALV.  ELECRAT* ; ELECRAT 
10B0  CN-TNSC4-T7nCHACR-T7\CAL*1/3 

1090  CHACP-<-<3*U-NA>/NA)t<CAL+2tT7>*7ST?)/<4-(3m-NA>/NA>> 

1199  CNCLI ■CN+CHACP 

1119  CHCL-CMCLI-T7 

1129  CALCL3a(NA/( 1-NA)  XCHCLI 

1159  IAL-1/3M7-T7MH-TH 

1169  IALCL3-NC* < <  T7+4/3 28NC- 1 ) > 

1179  BHCL- < 1-NC )*( (T7+4/3 >/<2*NC-l ) ) 

1189  BHAC*BMCL+BM 

1 199  ELECRATS • ( BALCL3* 1 33 . 34+BNCL* 146 . 64 )/ ( CNCLI * 1 46 . 64+CALCL38 133.34) 

1195  EXCLYRAT-ELECRAT-ELECRATS 

1299  EXCLY • EXCLYR AT 8 ( CNCL l t . 14664+CALCL3* . 13334 ) 

1229  PRINT\PRINT\PRINT 

1239  PRINT  'CIUEN  THE  ABOUE  PARAMETERS: ‘SPRINT 

1249  SPEN-26.8SU/(9.14664«<»nCL+CHCLI)«9.13334UBALCL3+CALCL3>+9.92698t(CAL>> 

1243  IF  EXCLV<9  THEN  SPENXCLV*0\GOTO  1259 
1245  SPCNXCLV-26 . 88U/ ( ( 26 . BtU/SPEN )+EXCLV ) 

1259  PRINT  'BALANCED  ANODE  REACTION  IS ‘SPRINT 

1269  PRINT  USING  '(tt.lUt  NCL  ♦  88.8888  ALCL3 )  ♦  8.999  AL  ♦  8.8888  ALCL4-  ♦  8.88 
88  NALCL4',  CNCL,CALCL3.CAL,C4,CNACR 

1265  PRINTS  PRINT  USING  ' - >  88.8888  NALCL4  ♦  8.8888  H*  ♦  1.9999  E-'.CHACP.CM 

1279  PRINTsPRINTsPRINT  'BALANCED  CATHODE  REACTION  IS'sPRINT 

1289  PRINT  USING  '(88.8888  NCL  ♦  88.8888  ALCL3 >  ♦  8.8888  N*  ♦  1.9999  E-',BNCL,BAL 
CL3.BN 

1299  PRINTSPRINT  USING  ' - >  8.888  AL  ♦  88.8888  HALCL4  ♦  8 . 8888AL2CL7- ' , BAL , BNAC, 

B7 

1292  PRINTSPRINTSPRINT  USING  'THE  STOICHIOMETRIC  RATIO  OF  HCLV/HALV,  ELECRATS.  IS 
88. 888 '.ELECRATS 

1294  PRINT  SPRINT  USING  'THE  SPECIFIC  ENERGY  IS  888.888  UHR/KG'.SPEN 

1296  IF  EXCLY<0  THEN  PRINTsPRINT  USING  'THE  STOICHIOMETRIC  RATIO  EXCEEDS  THE  SPEC 

IFIED  HCLV/MALY  RATIO,  8.88', ELECRAT 

1399  IF  EXCLV<0  GOTO  1392 

1391  PRINTsPRINT  USING  'SPEC  ENERGY  INCLUDING  EXCESS  CATHOLYTE  IS  98.89  UHR/KG'.S 
PENXCLY 

1392  PRINT  82,'2*,'iNAj',*;NCj',*;T7,', ‘ jTM; *» ')Wj  *« * jSPEN; • , • jSPENXCLY 
130S  NEXT  TH 

1307  CLOSE  82 

1310  END  ' 


Ready 


Fig  4  SPENCALC  program  for  calculating  specific  energy  for  i  concentration 
cell  with  fixed  transport  numbers. 


spencalc 


20-0CT-1983  14:54 


TRANSPORT  NUMBER  FOR  AL2CL7-,T7?  .25 
TRANSPORT  NUMBER  FOR  IMIDAZOtIUM  ION. 
MOLE  FRACTION  OF  ANOLVTE,  NA?  .37 
MOLE  FRACTION  OF  CATHOLYTE.  NC?  .60 
BATTERY  LOAD  UOLTAGE,  U?  1 . 1 
EXPERIMENTAL  RATIO,  MASS  CLV/MASS  ALY, 


TM7  .75 


ELECRAT?  1. 
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GIUEN  THE  ABOUE  PARAMETERS: 

BALANCED  ANODE  REACTION  IS 

(  3.5365  MCL  ♦  2.2532  ALCL3 )  ♦  0.333  AL  ♦  0.2500  ALCL4-  ♦  0.2500  MALCL4 
- >  3.0865  MALCL4  ♦  0.7500  M+  ♦  1.0000  E- 

BALANCED  CATHODE  REACTION  IS 

(  3.1667  MCL  ♦  4.7500  ALCL3 )  ♦  0.7500  M*  ♦  1.0000  E- 

- >  0.333  AL  ♦  3.9167  MALCL4  +  0.25OOAL2CL7- 


THE  SPECIFIC  ENERGY  IS  14.9664  UHR/KG  .. 

THE  STOICHIOMETRIC  RATIO,  MASS  CLV/MASS  ALY,  Is  1.272 

SPECIFIC  ENERGY  INCLUDING  EXCESS  CATHOLYTE  IS  13.6067  UHR/KG 
Ready 


Fig  5  Output  from  run  of  SPENCALC 


Table  1. 

Output  of  SPENPHD  calculation  with  parameters  fixed  to  match  those  of 
experimental  cells. 


FILE*  CASE37S60 
CASE37S60 


NA 

NC 

T7 

TM 

U 

SPEN 

SPENXCLV 

C3 

C3 

C3 

CD 

CU0LTS3 

CUHR/KGD 

CUHR/KGD 

0.37 

0.60 

0.50 

0.50 

1.1 

12.93 

11.76 

0.37 

0.60 

0.45 

0.55 

1 . 1 

13.29 

12.09 

0.37 

0.60 

0.40 

0.60 

1.1 

13.68 

12.43 

0.37 

0.60 

0.35 

0.65 

1.1 

14.08 

12.80 

0.37 

0.60 

0.30 

0.70 

1 . 1 

14.51 

13.19 

0.37 

0.60 

0.25 

0.75 

1 . 1 

14.9? 

13.61 

0.37 

0.60 

0.20 

0.30 

1.1 

15.45 

14.05 

0.37 

0.60 

0.15 

0.35 

1 . 1 

IS. 97 

14.52 

0.37 

0.60 

0.10 

0.90 

1 . 1 

16.52 

15.02 

0.37 

0.60 

0.05 

0.95 

1.1 

17.12 

15.56 

0.37 

0.60 

0.00 

1.00 

1 . 1 

17.76 

16.15 
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of  SPENXCLY  is  in  effect  simply  1.5  *  MALY,  and  since  MALY  is  constant  for  a 
given  NA,  SPENXCLY  is  also  constant.  (Obviously,  MALY  would  change  if  the 
number  of  equivalents  of  electrons  to  be  generated  were  changed  from  the  value 
of  1.0  used  in  these  calculations.)  Note  that  this  value  of  SPENXCLY  does  not 
apply  to  the  cells  which  have  NC  of  0.59  or  less.  These  runs  have  SPENXCLY' S 
set  at  zero  since  the  calculated  values  of  ELECRATS  for  these  cases  exceed  the 
specified  ELECRAT  of  1.5. 

Table  1  is  a  listing  of  the  results  of  a  SPENPHD  calculation  with  NA,  NC, 
V,  and  ELECRAT  values  which  were  used  in  actual  concentration  cells.  These 
data  are  plotted  in  Fig  8  along  with  the  experimental  data  reported  in  (2). 


IV.  Discussion: 


The  highest  theoretical  specific  energy  that  one  can  associate  with  the 
concentration  cell  discussed  here  would  be  that  obtained  with  pure  AlCl^  as 
the  catholyte  and  pure  MCI  as  the  anolyte.  In  such  a  case,  we  could  write  the 
following  cell  reactions: 


ANODE:  A1  +  4C1  ♦  A1C1.  +  3e 

4 


3e”  +  A1C13  A1  +  3Cl" 


CATHODE : 


3C1~+  3A1C1,  +  3  A1C1. ” 

_ _ 3  4 

3e”  +  4A1C13  ♦  3A1C14“  +  A1 

4  Cl"  +  4A1C13  +  4A1C14” 


Even  for  this  cell  the  specific  energy  is  only  98  Whr/kg  assuming  a  load 
voltage  of  1.4  V  as  in  the  cells  discussed  above.  Thus  the  cell  is  clearly 


♦  ‘  e  *  •  **0 


not  being  studied  as  a  potential  high  energy  density  battery.  These 
calculations  and  their  comparison  to  experimental  results  were  primarily 
intended  to  increase  understanding  of  the  electrochemistry  of  the  MC1/A1C1^ 
melts. 

The  experimental  results  in  Fig  8  are  in  reasonable  agreement 
with  the  calculated  results.  It  should  be  recalled  that  the  t+  values  of 
the  abscissa  are  true  transport  numbers  for  the  experimental  results  only. 
For  the  calculated  results  the  abscissa  values  are  actually  t+  +  t 
Thus  a  strict  comparison  of  the  calculated  with  experimental  results  would 
require  the  calculated  curves  to  be  displaced  to  the  left  by  an  amount  equal 
to  the  value  of  t_^.  However,  experimental  results  show  that  the  value  of 
t  ^  in  these  cells  is  probably  small,  a  value  of  -0.05  ±  0.18  having  been 
reported  (2). 


Ue  see  lower  SPEN  values  for  the  experimental  data  because  the 
discharges  of  actual  batteries  can't  be  carried  to  neutral  melt  formation. 

The  average  final  NA  in  the  experimental  cells  (2)  was  0.48,  whereas  in 
calculations  done  here,  the  final  NA  is  0.50  (with  NA  *  0.37  initially  in  both 
cases).  Thus  there  is  strong  evidence  that  the  model  of  the  concentration 
cell  used  in  these  calculations  is  fundamentally  correct. 


V.  Conclusions: 

We  can  calculate  the  specific  energy  for  MC1/A1C1^  concentration  cells  as  a 
function  of  electrolyte  concentrations  and  transport  numbers.  The  calculated 
values  are  in  good  agreement  with  previously  obtained  experimental  data, 
supporting  the  model  on  which  the  calculations  are  based. 
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’ig  8  Plotted  results  of  SPENPHD  calculation  for  cell  with  parameters  fixed 
to  match  experimental  cells  compared  to  experimental  results,^ 
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